Detecting planets in protoplanetary disks: A prospective study 
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ABSTRACT 

We investigate the possibility to find evidence for planets in circumstellar disks by infrared and 
submillimeter interferometry. We present simulations of a circumstellar disk around a solar-type 
star with an embedded planet of 1 Jupiter mass. The three-dimensional (3D) density structure 
of the disk results from hydrodynamical simulations. On the basis of 3D radiative transfer 
simulations, images of this system were calculated. The intensity maps provide the basis for the 
simulation of the interferometers VLTI (equipped with the mid-infrared instrument MIDI) and 
ALMA. While MIDI/ VLTI will not provide the possibility to distinguish between disks with or 
without a gap on the basis of visibility measurements, ALMA will provide the necessary basis 
for a direct gap detection. 

Subject headings: hydrodynamics, radiative transfer — techniques: interferometric — (stars:) circum- 
stellar matter, planetary systems, pre-main sequence 

1. Introduction 

Based on studies of the evolution of protoplan- 
ets in protoplanetary disks, it has been established 
that - depending on the hydrodynamic properties 
of the planet and the disk - giant protoplanets 
may open a gap and cause spiral density waves 
in the disk (see, e.g., Kley 1999, Kley et al. 2001). 
The gap may extend up to several AU in width. 



Thus, the question arises if one can find such a 
gap as an indicator for the presence of a proto- 
planet with present-day or near-future techniques. 
In order to study this possibility, we use hydrody- 
namical simulations of a protoplanetary disk with 
an embedded planet (§ 2) and compute the ex- 
pected brightness distributions, with a 3D radia- 
tive transfer code (§3). The resulting images are 
presented in § 4. Finally, we show that ALMA 



will provide the necessary basis to detect gaps in 
circumstellar disks in the mm/submm wavelength 
range (§ 5). 

2. The disk model 

The density structure of the protoplanetary 
disk - on which the subsequent radiative trans- 
fer (RT) simulations and discussions are based 
on - results from hydrodynamical simulations in 
which the disk is assumed to be flat and non- 
sclf-gravitating. The mutual gravitational inter- 
action between the planet and the central star, 
and the gravitational torques of the disk acting 
on planet and star are included. The 3D den- 
sity structure is Gaussian in the vertical direction 
where for the scale height H(r) we assume a con- 
stant ratio H/r=0.05. The radial density profile of 
the disk is shown in Fig. 1[A]. 
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Fig. 1 . — Radial density profile in the midplane of 
the disk (averaged in <j> direction, arbitrary scale). 
The thin lines represent the corresponding density 
distributions without an embedded planet. [A] 
With accretion onto the central object (open inner 
disk boundary). [B] No accretion onto the central 
object (closed inner disk boundary). 

The mass of the star is assumed to be 1 M , the 
mass of the planet is 1 Jupiter mass, and the mass 
of the disk is 0.05 M Q . For the dust-to-gas mass 
ratio we take the classical value of 1:100. The di- 
ameter of the disk is 104 AU which is comparable 
to the size of the disks seen in silhouette against 
the Orion Nebula (McCaughrean & O'Dell 1996). 
The planet is located at a distance of 5.2 AU from 
the star. The structure of the spirals and the 
gap reaches an equilibrium after about 150 orbits 
(« 1800 yrs). The specifics of the models and the 
results, in particular the width and internal struc- 
ture of the gap, are described in detail in Klcy 
(1999) and Kley et al. (2001). 



3. Radiative transfer simulations 

We used the results from the hydrodynamical 
simulations to calculate the dust continuum emis- 
sion. The RT has been performed with a 3D con- 
tinuum code which is based on the Monte Carlo 
method. The RT is simulated self-consistently, 
taking into account both the initial temperature of 
the dust due to viscous heating and the additional 
energy input of the central star. The dust density 
and temperature structure are defined on a grid 
which is chosen to be identical with that of the hy- 
drodynamical simulations. Thus, any additional 
discretization error resulting from the "transfer" 
of the density/temperature structure to the RT 
code is avoided. For further information concern- 
ing the RT code, we refer to Wolf et al. (1999) and 
Wolf & Hcnning (2000). 

In addition to the results from the hydrodynam- 
ical simulations, we introduce the following RT pa- 
rameters: spherical dust grains consisting of "as- 
tronomical" silicates (optical data from Draine & 
Lee 1984, radius 0.12/im); star: effective temper- 
ature T e ff = 5500 K, luminosity L = 1L©; wave- 
length range for the simulation of the radiative 
transfer: 0.03 . . . 2000 fim. We should note that 
the resulting temperature structure of the disk is 
in good agreement with the temperature distribu- 
tion resulting from 2D hydrodynamical disk mod- 
els with the same geometry. 

4. Images 

In Figure 2, the images of the inner region (di- 
ameter 28 AU) of the disk, seen under different in- 
clination angles i (i = 0° - face-on, 30°, and 60°), 
at A=700^m are shown. Assuming a distance of 
140 pc of the object (which is the distance to the 
nearby star- forming region in Taurus), the pixel 
scale is 2 mas/pixel. The gap can be clearly seen 
- even in the case of large inclinations of the disk. 
The spiral density waves - which are also caused 
by the motion of the planet - are difficult to trace. 

Due to the huge intensity peak caused by 
warm/hot dust located very close to the star 
(1. . .2 AU), the outer regions of the disk remain 
dark at near/mid-infrared wavelengths. The gap 
located at 5AU can be hardly seen. In contrast 
to this, in the sub- millimetre/millimetre wave- 
length range the brightness distribution is much 
smoother. This can be explained by the be- 
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Fig. 2. — Images of the inner region (diameter 28 AU; see Fig. 1 1[A] for the corresponding radial density 
profile) of the circumstellar disk with an embedded planet of 1 Jupiter mass at a wavelength of A = 700 /im 
and inclinations i = 0° (face-on), 30°, and 60° (see symbols in the bottom row). Each frame consists of 

101 x 101 pixels. The gap is clearly visible since the density drops by two orders of magnitudes in this region 

i 

(optical depth at a wavelength of 700 /im: 2 x 10~ 3 ). To increase the color dynamic range, the quantity 1^ 
is shown. The spiral density waves in the disks are weakly visible at i — 0°. The median specific intensity 
h [W/im^sr- 1 !!!- 2 ] per pixel amounts to 1.74-10 - 23 (i=0°), 1.68-10 - 23 (i=30°), 9.93-10~ 24 (i=60°). See 
electronic edition for colored images. 



haviour of the Planck function B\(T). The ratio 
B x {Tx) / B\{T 2 ) (Ti > T 2 ) strongly decreases with 
increasing wavelength. Thus, the flux ratio and 
the brightness contrast between regions character- 
ized by different dust temperatures (Ti and T 2 ) 
decreases with increasing wavelength A. 

5. Simulations of observations with MIDI 
(VLTI) and ALMA 

In Figure 3 visibility curves at A=10/mi for 
two different inclinations of the disk are shown 
(with/without gap). The goal of these simula- 
tions is to check whether MIDI at the VLTI can 
be used to detect gaps. In order to avoid sam- 
pling problems in the image plane which could 
affect the calculation of the visibilities, the res- 
olution of the images had to be increased by a 
factor of 4 along each coordinate axis (compared 
to Figure 2). It was found that the visibilities at 
a given baseline differ by less than 5%. This is 
due to the fact that the density distribution in 
the innermost region - which dominates the 10 /im 
flux - is only negligibly affected by the presence of 
the planet at a distance of 5.2 AU. Taking into ac- 
count uncertainties of "real" measurements (e.g. 
with MIDI - the mid-infrared intcrfcrometric in- 
strument for the VLTI), a distinction between dif- 



ferent disk models (with/without gap) based on 
significant differences between the visibility pro- 
files is not possible even by a beam combination 
of the most distant telescopes of the VLTI. To 
demonstrate the influence of the particular den- 
sity profile of the disk, in Fig. 4 the visibilities of 
a disk, with the same parameters as described in 
§ 2 but without accretion onto the central star, are 
shown (see Fig. 1 [B] for the radial density profile) . 
While again the presence of a gap has almost no 
influence on the visibility, the visibility profile is 
flatter compared to the visibilities shown in Fig. 3. 
Thus, one could clearly distinguish a disk with ac- 
cretion from a disk without. 

In contrast to a two-beam interferometer such 
as MIDI at the VLTI, ALMA (Atacama Large Mil- 
limeter Array) will combine signals of 64 anten- 
nas with an aspired maximum baseline of 12 to 
14 km. It will cover the submm/mm wavelength 
range (v—iQ. . .900 GHz). Due to the large num- 
ber of receivers and the resulting broad distribu- 
tion of baselines, a sufficient u-v plane coverage 
can be achieved after a few hours of observation 
and image reconstruction will be possible. In addi- 
tion, and in contrast to the first instruments at the 
VLTI, the visibility phase can be easily measured 
in the millimetre wavelength range. In Figure 5 a 
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Fig. 3. — Normalized visibilities at A=10/im for a 
disk with an inclination of i — 0° and 60° (right 
column; assumed distance: 140 pc). For compar- 
ison the same profiles in the u-v plane but for a 
disk without a gap (but with the same mass and 
inclinations) are also given (left column) . The vis- 
ibilities marked by the thick and the thin line are 
oriented perpendicular to each other in the u-v 
plane (parallel to the major/minor axis of the el- 
lipse resulting from the projection of the disk onto 
the plane of the sky). 



reconstructed image of the disk based on a simu- 
lation of ALMA observations is shown (disk seen 
face-on, A=700/im). Even under consideration of 
the thermal noise caused by a system temperature 
of T sys =500K, the gap is clearly visible. 

6. Conclusions 

Based on hydrodynamical simulations and sub- 
sequent 3D continuum RT calculations we gener- 
ated millimetre images of a circumstellar disk with 
an embedded planet surrounding a solar-type star. 
The gap can be seen very clearly in the simu- 
lated images but an extremely high angular res- 
olution is required (« 10 mas). Because of the ex- 
treme brightness contrast in the innermost region 
of the disk in the near to mid-infrared, the gap can 
hardly be detected in this wavelength range with 
imaging/interferometric observations. 

We found that it will be possible to distinguish 
between a disk with or without accretion onto 
the central star with MIDI, but the visibility at 
10 /urn is almost insensitive to the presence of a 
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Fig. 4. — Normalized visibilities at A=10/mi for a 
disk with an inclination of i = 0° and 60° (right 
column; assumed distance: 140 pc). In contrast 
to the visibilities shown in Fig. 3, a disk without 
accretion onto the central star is assumed here. 



gap. In contrast to this, the millimetre interfer- 
ometer ALMA which will become available in the 
near future will achieve this goal: it will provide 
the basis for the reconstruction of an image of a 
gap. Thus, the search for massive protoplanets in 
circumstellar disks can be based on the indication 
of a gap. For this purpose, the longest possible 
baselines are required. 

This research was supported by the DFG grant 
Kl 650/1-1. We appreciate the suggestions given 
by G. Brydcn who was the referee of this article. 
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Fig. 5. — Reconstructed image of the disk seen 
face-on resulting from a simulation of ALMA. The 
gap at an angular distance of 37 mas (5.2 AU) from 
the star is clearly visible. Wavelength: 700 jitm 
(428 GHz); bandwidth: 8 GHz; total integration 
time: 4 h; system temperature: 500 K; phase noise: 
30°; max. baseline: 10km. 



This 2-column preprint was prepared with the AAS IATjrjX 
macros v5.0. 
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